SU-4240-691 



The Renormalization Group and its Finite 
Lattice Approximations 



Angelo Cacciuto * Eric GregoryWd Alex Travesset* 

Physics Department, Syracuse University, 
Syracuse, NY 13244-1130, USA 



Abstract 

We investigate finite lattice approximations to the Wilson Renor- 
malization Group in models of unconstrained spins. We discuss first 
the properties of the Renormalization Group Transformation (RGT) 
that control the accuracy of this type of approximations and explain 
different methods and techniques to practically identify them. We also 
discuss how to determine the anomalous dimension of the field. We 
apply our considerations to a linear sigma model in two dimensions in 
the domain of attraction of the Ising Fixed Point using a Bell- Wilson 
RGT. We are able to identify optimal RGTs which allow accurate 
computations of quantities such as critical exponents, fixed point cou- 
plings and eigenvectors with modest statistics. We finally discuss the 
advantages and limitations of this type of approach. 
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1 Introduction 



1.1 Scope and organization of the paper 

Finite lattice approximations provide an alternative and, we think, largely 
unexplored approach to the Renormalization Group (RG) [1]. The advantage 
of a direct RG approach lies in its generality as not only does it provide, in 
principle, a direct and complete calculation of the critical exponents of the 
model, but also further valuable information such as the fixed point Hamil- 
tonian and the renormalized trajectory. In practical calculations, however, 
the apparent simplicity of the method is usually complicated because one 
needs to introduce further approximations (truncation in the operator basis, 
too few iterations of the RG, etc.). In addition there are errors inherent to 
the approach (finite size effects, statistical errors in numerical simulations, 
etc.). Models of unconstrained spins present, in addition, the problem of 
determining the rescaling of the field. 

These difficulties may be overcome if a detailed knowledge of the prop- 
erties of the renormalization group transformation used are available. We 
discuss these properties in some detail and explain different criteria, some 
already known in the literature, to single out the best renormalization group 
transformation and determine the anomalous dimension of the field. 

This paper is organized as follows. In Sect. 2 we discuss first some general 
aspects of the Renormalization Group (RG) and the Monte Carlo Renormal- 
ization Group (MCRG) technique. We then describe the parameterizations 
of the effective Hamiltonians generated along the RG flow more appropriate 
for the kind of approximation to be applied. We finally introduce the partic- 
ular RGT to be used in this paper and we discuss the physical meaning of 
the parameters involved. 

In Sect. 3 we address the problems and subtleties that appear in the finite 
lattice approximations of the RG. We first discuss the theoretical aspects and 
then explain different methods and techniques to be used in practice. We 
also discuss on the different sources of possible errors involved. 

The technical details of the algorithms used, as well as the description of 
the actual parameters of the simulation, are presented in Sect. 4. 

In Sect. 5 we present in detail the numerical results concerning short- 
rangeness of the effective interactions generated along the flow, optimization 
of the transformation, determination of the anomalous dimension, calculation 
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of the critical exponents and the estimate of the different errors involved in 
the calculation. 

The results presented in Sect. 5 are rather technical, so we summarize 
them in the first part of the conclusions. In the second part we discuss further 
improvements and future outlook. Some technical issues are discussed in the 
appendices. 

2 The RG 

2.1 Some generalities of the RG 

The idea of the RG[1] is simple; average out short-distance modes leaving 
large distance properties intact. Let us assume that our model of interest 
is described by the Hamiltonian (so called canonical surface in the RG lan- 
guage) 

w°M = E*2o«(0> (i) 

a 

where {Kp}p =1 ... are the couplings and {O a }a=i, - an arbitrarily chosen basis 
of operators. The integration of short- distance modes is implemented by a 
transformation 1Z, the Renormalization Group Transformation (RGT), whose 
result is to define a new effective interaction Ti 1 which we parameterize with 
the same basis of operators as the canonical surface. The net effect of the 
RGT is a transformation in the space of all possible couplings of the model, 

K = n a {{KD) . ( 2 ) 

Iteratively applying Eq. 2 we generate a sequence of points in coupling space 
that define a trajectory. We can visualize, then, the effect of a particular 
RGT as trajectories flowing in coupling space. Usually, these trajectories 
converge to a special point, say {K*} a= i t ..., a fixed point(FP). It is a fixed 
point of the RGT Eq. 2, 

K* a = n a ({K;}) . (3) 

The importance of a FP speaks for itself. The universal properties of all 
Hamiltonians lying in trajectories converging to this FP are completely char- 
acterized by the the RG flow in an infinitesimal vicinity around it [1] ; all the 
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critical exponents, for example, just follow from diagonalizing the matrix 



-a/3 



en a (K) 



(4) 



{K*} 



In this way, the RG provides a completely general and exact prescription 
to compute the critical properties of any model. 



2.2 The MCRG method 

The exact T matrix is a difficult object to compute. There are a number 
of approaches available, but as our interest lies in applying numerical tech- 
niques, we seek an expression for the T matrix as a function of expectation 
values of operators. Remarkably, this can be achieved [2]. At the (i + l)th 
iteration of the RGT Eq. 2 one can easily prove [2] 



dll a {K) 



with 7 running over all possible operators and 



= £H^M 7 /3> (5) 



M 7 , = (0^0f)-(0^)(0f), 

H 7/3 = (0^0^)-(0^)(0^) . (6) 

These are expectation values of operators at the ith and (i + l)th iteration 
of the RGT, but can be computed using H k with k < i. If W is the FP- 
Hamiltonian, which obviously is the case if % = oo, Eq. 5 is the exact T 
matrix, as defined in Eq. 4. 

With the very same effort, the couplings of the effective interaction 7i^ +1 ^ 
can also be computed [3]. The formula is 1 

^ = -r£ H ^XA>, (7) 

a P a 

where the H matrix is defined in Eq. 6, A a = (0^ +1 ), and d a is the dimension 
of the operator O a , to be defined precisely later. 

1 This formula applies for linear sigma models. Other models require some adaptations. 
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2.3 The effective interactions: Organization of the ex- 
pansion 

We are interested in models of unconstrained spins. The simplest case of 
this type of model is the linear sigma model. This is the model to which we 
explicitly apply all our considerations. The canonical surface is 

^°[0] = E{-^HE^ + ^) + ^) 2 + M0(n) 2 -l) 2 } , (8) 

n { fi ) 

where n runs all over sites of the lattice, which we take as an infinite simple 
square lattice, and \i runs just over nearest neighbors. Spins are uncon- 
strained, i. e. <f>(n) G (—00,00). 

Under iteration of RGTs, the new effective interactions 7i l will certainly 
contain more operators than those three in Eq. 8. In general, after i-RGTs 
we have 

-W[4\ = Y,KOa{*) , (9) 

a 

with a running over all possible operators compatible with the symmetries 
of the model. We need an efficient organization of the operators appearing 
in Eq. 9. 

We consider a basis of operators constructed out of monomials centered 
at an arbitrary site n, consisting of powers of fields at site n and products 
with fields in other lattice sites. A typical monomial is 

M k , K {n) = E 0H 2fe 0(n + £)••• 0(n + w) , (10) 
{u,---,w}e/c 

where /C is the class of all vectors that starting from some representative ones 
may be constructed out of the symmetries of the lattice. The full operator 
is constructed by summing over all lattice sites, 

O k ^)=T. M Kdn) ■ (11) 

n 

We classify all possible monomials (and therefore the operators constructed 
out of them), according to their type, length (l a ) and dimension (d a ); 

• The type: We define 4 different types, shown explicitly in Table I. Type 
are bilinear operators. Type 1 are even powers of the lattice field at 
point n. Type 2 and 3 include products of bilinears and even powers 
of the field at site n. 
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Type u d a I, 



<f>(n)(f>(n + u) 


(1,0) 


2 


1 




(1,1) 


2 


1 




(2,0) 


2 


2 




(2,1) 


2 


2 




(2,2) 


2 


2 




(p,q) 


2 


P 


Type 1 


<f)(n) 2k 




2k 





Type 2 


u 






<f>{n)<f>{n + u)<f)(n) 2k 


(p,q) 


2k+2 


V 


Type 3 


u, w 






4>{n + w)(f)(n + u)(j)(n) 2k 


(p,q) , (r,s) 


2k+2 





Table 1: Operators included in the expansion, in the order they are consid- 
ered with their associated dimension and length. The vectors u = (p, q) are 
the representative vectors of the class /C in Eq. 10. The parameter k starts 
at 1. 

• The length l a : This is the maximum linear separation between fields 
in a given monomial. An operator with length l a can be accomodated 
in a lattice as small as (2l a + l) d , d being the dimensionality of the 
system. 

• The dimension d a : This is the number of product fields in each mono- 
mial. 

We organize operators within a certain type in increasing length and 
dimension. Overall, we are performing a systematic expansion with the 
assumption that the RG-couplings decay rapidly both with dimension and 
length. We elaborate more on this in considering the properties of RGTs. 

There is a last remark concerning the canonical surface Eq. 8. In this 
paper, the model will be considered in 2 dimensions. The most general RG 
flow is rather complicated as there are an infinite number of inequivalent FPs. 
From the canonical surface Eq. 8, however, we just can reach two of them, 
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the Gaussian FP (GFP), which is the most infrared unstable, and the Ising 
FP (IFP), which is the most infrared stable. To reach any other FP, higher 
operators, together with a fine tuning of the couplings, should be considered. 

2.4 The Bell-Wilson RGT 

There are (infinitely) many different RGTs giving rise to different RGTs 
having different properties. The utility of those properties lies in that they 
may be used to single out some RGTs as being more accurate than others in 
approximate calculations of the T matrix. The choice of a particular RGT 
is therefore determined by the type of approximation we apply. 

As we are eventually interested in a finite lattice approximation, a very 
important property is to identify RGTs generating very short-ranged effective 
Hamiltonians. That is, at each RG iteration, the effective interaction H. 1 
generated can be parametrized with a few operators of small dimensions and 
lengths. This is the case if FP couplings, for example, decrease at least as 
rapidly as 

\K*\<Ce- da ^ e -i<*/ti , (12) 

i./e./ exponentially with both dimension and length. We introduced two 
decay corrrelation lengths £z and £d- The importance of this property lies 
in that the more short-ranged the FP, the more insensitive it will be to 
truncation of higher dimension and length operators. This is the rationale 
behind the organization of the expansion of the operators generated along 
the RG flow in the previous section. 

The RGT first introduced in [4] has the property that the Gaussian FP 
is short-ranged, in the sense of Eq. 12 [4] [5] [6]. It is therefore a natural 
candidate to investigate for other FPs as well. The RGT is 

e -n^[#] = I Y[dcj>(n) e- n ^ WnflME ^B ^ ? (13) 

J n 

where the average is over c x c cells in the fine lattice, with c = 2 as shown 
in fig 1. There are two a priori free parameters, aw and b. 

The parameter aw just labels different transformations. Universal quan- 
tities are independent of it. Of course, within an approximation that may 
not be the case exactly, and it provides us with a freedom that we can use 
to optimize the performance of the RGT. 
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n 







c= 2 

Figure 1: Graphical representation of the short-distance integration of the 
RGT Eq. 13, n labels sites in the fine lattice, % in the coarse one. 

The parameter b, however, has a much deeper meaning related to the 
unboundeness of the spins; we may rescale them, and yet the physics of the 
model is the same. If there is a FP we may generate a whole line of new 
FPs by rescaling the field 0. The model shows a redundant direction [1], 
[10]. By definition, any point in that line is a FP of the transformation, 
so the RGT cannot make the FP flow along this direction, and therefore it 
is exactly marginal. A first consequence of this is that the T matrix will 
always have at least one eigenvalue exactly equal to one. The parameter b 
gets uniquely and self-consistently fixed by requiring that there is no flow 
along this marginal redundant direction. Any other choice of b would cause 
this line to be relevant or irrelevant, and there would be no FP in the theory, 
in contradiction with the physics of the model. The rescaling parameter is 
not a freedom of the model. In fact, it has a universal meaning. It may be 
related to the anomalous conformal dimension of the field rj, the exponent 
controlling the universal algebraic decay as a function of distance of the two- 
point correlation function 



G{x) = (0(^)0(0)) 



\rj,\d—2+rj 



X 



OO 



(14) 



with 



rj = d + 2 -2 



In b* 
In c 



(15) 



where b* stands for the rescaling factor b at the FP and c will be usually 
equal to 2 [1]. 
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3 The RG in a finite lattice 



3.1 Theoretical issues 

We consider now the canonical surface Eq. 8 and the transformation Eq. 13 
defined on a finite lattice with periodic boundary conditions (so that trans- 
lation symmetry is preserved). The T- matrix may be computed from Eq. 5 
and Eq. 6. The question is how much the T-matrix computed within this 
approximation differs from its infinite volume value. 

Let us assume that the FP of the RGT is such that the largest operator 
length is n s . The entire FP may be accommodated in a lattice as small as 
(2n s + l) d . Let us consider it now in a (2(2n s + l)) d volume, and apply 
a RGT, reducing the volume to a (2n s + l) d . As the FP still perfectly 
fits in, the FP couplings should not differ from their infinite volume values. 
There are no finite size effects 2 . This is even more remarkable when one 
considers that both the couplings and the T-matrix are computed from the 
matrices Eq. 6. These involve the computation of expectation values that 
show severe finite size effects as the system is critical, and hence has infinite 
correlation length. This may be one of the most interesting features of MCRG 
methods. We may, at least in principle, obtain infinite volume results working 
within a finite lattice approximation. Finite size effects are not dictated by 
the correlation length (which is infinite) but by the range of the FP (and 
therefore of the RGT). A consequence of the previous arguments is that the 
parameters of the canonical surface must be fine tuned at criticality. We 
need to perform the calculations at the values of the exact critical surface of 
the infinite dimensional system, at virtually infinite correlation length. 

Obviously, applying a RGT in a finite lattice reduces its volume. Even if 
the RGT is very short-ranged, it may happen that the FP cannot be reached 
before the lattice is so small that huge finite-size effects, as explained, enter 
into the game. As we may iterate the RGT only a finite number of times, 
we need RGTs such that the canonical surface is very close to the FP [7] [8]. 

In short, any RGT to be used in a finite lattice approximation should 
fullfill two conditions; 

• Generate short-ranged effective Hamiltonians. 

2 More Rigorously, finite size effects are exponentially suppressed and therefore very 
small. They are irrelevant for this argument. 
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• Approach the FP rapidly. 

In models of unbounded spins there is still the issue of the determination 
of the rescaling factor b in the RGT Eq. 13. This parameter is intimately 
related to the existence of a line of FPs. In general, however, not all FPs 
along this line will be short-ranged enough and the existence of the whole 
line may be apparent only in going to prohibitively large lattices. However, if 
a sufficiently short-ranged FP exists, we may expect at least an infinitesimal 
portion of this line to survive around it and the T matrix to exhibit at least 
one eigenvalue equal to one with good precision. 

3.2 RG in practice 

The critical couplings of the canonical surface Eq. 8 in an infinite lattice can- 
not be computed within a MCRG approach. It is a quantity to be computed 
from other techniques, and we extracted it from the analysis of existing coef- 
ficients of the strong coupling expansion of the canonical surface, as discussed 
in appendix A. 

Once the critical couplings are known, iteration of the RGT Eq. 13 would 
drive the canonical surface towards the most infrared stable FP, the Ising 
FP (IFP) in our model. The strategy we follow is to pick a value for a w 
(see Eq. 13) and apply as many RGTs as allowed (that is, until the lattice 
becomes a 4 2 volume) with different values for the rescaling value b (as defined 
in Eq. 13). In the previous subsection we discussed the different issues that 
need to be addressed. The different criteria we used to tackle them are as 
follows, 

1. Short-ranged effective Hamiltonians: 

(a) Using Schwinger- Dyson Eq. 7 we compute the couplings and check 
for the ansatz Eq. 12. 

2. Rapid approach to the FP: There are two different criteria to be met. 

(a) If a FP is reached, the couplings of these effective Hamiltonians 
after successive applications of RGTs should coincide. 

(b) We perform RGTs starting from the same canonical surface but in 
a smaller lattice and compare expectation values of operators at 
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the same volume size. If a FP is reached these expectation values 
should agree. 

3. The determination of the anomalous dimension rj 

(a) The T matrix should exhibit an exactly marginal eigenvalue at 
the right value of the rescaling factor. 

(b) Similarly as in criterion 2(b), the expectation values of operators 
should agree at the correct value of the anomalous dimension. 

We then scan over different values of aw and select the values optimally 
fulfilling all previous criteria. 

Criteria 1(a) and 2(a) are obvious. Criteria 2(b) uses the fact that al- 
though finite size effects for these expectation values are huge, they enter in 
exactly the same functional form [9], and if a FP is reached, those expecta- 
tion values should agreee. Criterion 3(a) to pick up the anomalous dimension 
is a modification of [4]. It was first introduced in [11] and successfully ap- 
plied to the linear sigma model in 3 dimensions, although the lattices used 
were slightly too small, and finite size effect errors were difficult to assess. 
Criterion 3(b) is a necessary consequence of the discussion in criterion 2(b) 
and the fact that a wrong choice for b should not lead to a FP at all. 

Besides the statistical errors inherent in any numerical simulation, there 
are different sources of putative systematic errors, 

• e(l) Truncation errors: If some operator having a sizeable coupling is 
not included in Eq. 5 and Eq. 7 the H and M matrices miss sizeable 
matrix elements which translates into systematic errors in the determi- 
nation of both the T matrix and the couplings. 

• e(2) Finite size effects: This error appears when the lattice is too small 
to accommodate operators having sizeable couplings. 

• e(3) Off-criticality: If the couplings in the starting Hamiltonian are not 
tuned to criticality, the flow eventually converges towards the renormal- 
ized trajectory, moving away from the FP. 

• e(4) Lack of eigenvector convergence: It may happen that some eigen- 
vectors of the T matrix are just not convergent, in the sense that the 
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coefficients denning it have just a formal meaning. This may translate 
into the correponding eigenvalue being completely unreliable. 



• e(5) Redundant operators: Some eigenvalues may depend on the trans- 
formation, the associated critical exponent is then just an artifact of 
the transformation [12], [13]. 

Finally, let us recall that in models of unconstrained spins, like the ones 
we are interested in, the rescaling itself will have an associated error bar due 
to statistical errors and e(l),e(2) and most importantly e(3). This rescal- 
ing error bar the will translate into an additional error in computing other 
quantities. 

4 Numerical details 
4.1 The algorithm 

We use an embedding algorithm [14], each step consisting of 15 cluster- Wolff 
updates [15] and 5 metropolis hits. We need 20 sweeps of this type to com- 
pletely decorrelate data. We computed the error as a function of bin-size 
using standard Jackknife techniques and found a flat plateau showing that 
data in different bins are totally uncorrelated. After thermalization, we gen- 
erated 64,000 decorrelated 64 2 lattice configurations and 50,000 34 2 lattice 
configurations at criticality, and stored them in disk. The uncompressed 
storage size of these is about 2GB and 500MB respectively. The simulation 
was carried out on a 233 MHz Pentium processor, and took about 2 days of 
total CPU time. 

The RGT Eq. 13 is implemented using the shift 



where ( is a randomly distributed gaussian variable with zero mean and 
variance 1. From Eq. 16 we generate the successive RGTs. On a 233 MHz 
Pentium processor it takes about 100 minutes of CPU time to generate the 
M matrices and H matrices including up to 40 operators for 32 2 , 16 2 , 8 2 and 
4 2 volumes. 




c 



(16) 
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4.2 The method 



The canonical surface Eq. 8 is a actually a line k(\). We performed our 
simulations at A = 1.0, with corresponding k(1.0) = 0.6795. 

Concerning the expansion of the effective Hamiltonians generated along 
the RG flow, we considered up to 40 operators in total. All operators of 
type up to length 3, operators of type 1 up to dimension 10, operators of 
type 2 up to length 3 and dimension 6, and operators of type 3 up to length 
2 and dimension 4, as defined in Table 1. In the discussion that follows 
we will identify these operators by number from 1 to 40. Operators 1 to 5 
correspond to type 1 k = 1, ...5. Operators 6 to 14 correspond to type 
operators (p,q)= (1, 0), .., (3, 3). Operators 15 to 23 correspond to type 2 
operators with k = 1 (p,q)= (1, 0), .., (3, 3). Operators 24 to 32 correspond 
to type 2 operators with k = 2 (p,q)= (1, 0), .., (3, 3). Operators 33 to 40 
are the eight distinct k = 1 type 3 operators that fit in a 3 x 3 section of 
lattice (n at the center): [(p,q)(r,s)] = [(0,1)(0,1)], [(0,1)(1,0)], [(0,1)(0,-1)], 
[(0,1)(1,1)], [(0,-l)(l,l)], [(1,1)(-1,1)], and [(1,1)(-1,-1)]. 

We systematically explored anomalous dimensions rj = 0.10, 0.20, 0.25, 
0.30, 0.40, 0.50. In some particular cases additional values were also con- 
sidered. To improve the performance of the transformation, we studied the 
following values of the RGT parameter: a w = 8, 16, 20, 40, 80, oo. 

In order to efficiently handle the enormous amount of data generated, 
we organized the information in different html-tables, each one displaying 
the particular property under study. Samples of tables are available at the 
address http:/ /web. syr.edu/~acacciut 

5 Results 

In this section we present an extensive and detailed analysis of the methods 
and techniques already introduced. We keep the discussion technical. The 
main results may appear a little dispersed, so in the next section we provide 
a presentation of the most relevant final results and the different implications 
they may have. 
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5.1 Short-ranged Hamiltonians 



The effective interactions generated along the RG flow depend on aw (the 
transformation) and b (the rescaling). The subject under study now is the 
decay of couplings with dimension and length as a function of aw for a 
given b. The optimal RGT will be the one generating the most short-ranged 
effective interaction. 

The first issue to address is the dependence of couplings on length. All the 
different values of aw examined conform to the decay ansatz Eq. 12, although 
the decay is too fast to allow a precise determination of the correlation length 
decay £j, as coupling constants are zero within statistical error bars for I > 2. 
The typical situation ocurring for any value of a w is plotted in 3 fig. 2. 
From the data in fig. 2 we can give a very conservative upper bound for the 
correlation length decay, 



valid for any value of aw- We conclude that operators have a very good 
decay property for any aw value as far as length is concerned. 

Couplings as a function of dimension show a very strong dependence on 
a w- In fig- 3 a logarithmic plot of the magnitude of type 1 couplings is shown 
as a function of dimension. For small aw the decay ansatz, Eq. 12, sets in for 
small values of dimension. For aw = 8, we extract a dimension correlation 
length decay 



The correlation length decay as a function of distance is slower than with 
length. Increasing aw further, there is a transient increase of couplings, 
with a peak moving forward in dimension as a function of aw, before the 
asymptotic form Eq. 12 becomes apparent, as shown in fig. 3. In particular, 
the RGT at aw = oo shows a remarkably poor behavior. 

We cross-check previous results by studying the decay of type 2 operators. 
As shown in fig. 4, couplings decay exponentially, with a correlation length 
completely compatible with the bound in Eq. 17. 

From all the evidence accumulated, it is clear that smaller values aw re- 
strict the dependency of the effective interactions to lower dimension field 

3 We generalize slightly the definition of length by giving length y/2 to the operator 
(1,1) instead of 1, as defined. 




(17) 



& ~ i.o . 



(18) 
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Figure 2: Plot for aw — 16, r\ — .30 of the magnitude of type (bilinear) 
couplings as a function of length for 3 RGTs. Squares are vertical or horizon- 
tal products and circles the diagonal ones. The solid line is an exponential 
decay. 

operators. In particular, effective interactions generated with RG transfor- 
mation with aw > 80 are not particularly recommended, as the couplings 
decay too slowly with dimension to obtain reliable results. 

5.2 Fast approach to the FP 

The effective interactions show a double dependence on the rescaling b and 
on the parameter aw The dependence on b is related to the determination of 
the anomalous dimension and will therefore be studied in the next subsection. 
Now we content ourselves with identifying the RGTs that converge rapidly 
to putative FPs. 

Concerning criterion 2(a), the matching between couplings in successive 
RGTs, although the amount of information that may be extracted is certainly 
limited by the somewhat large statistical errors in the couplings, it is still 
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5.0 




Figure 3: The logarithmic plot of the magnitude of type 1 couplings against 
dimension, d. Results are for 3rd RGT at rj = .30. 



powerful enough to provide strong evidence that we approached a FP after 
a few RGTs. 

For values of aw > 80 couplings from successive RGTs do not match 
satisfactorily, and we cannot identify a FP for any value of the rescaling. In 
agreement with previous results, RGTs with large aw values perform very 
poorly. We move on to investigate lower values for aw- Results at aw = 40 
do not show a clear FP for any value of the rescaling, as there are some 
couplings that do not agree within 2a (statistical), see table 2. 

For aw = 8, there is no clear evidence for a FP at rj > .30, see table 3, 
although we find an acceptable matching for smaller values of rj. 

The value of aw = 16 shows putative FPs within the statistical errors 
for values of the anomalous dimension in the range ?7 = .20 — .40. As an 
example, in fig. 5 couplings from 2 and 3 RGTs are plotted at rj = .25. It 
is apparent that most of the couplings agree very well. It seems then that 
aw = 16 is an optimal choice according to criterion 2(a). 

Concerning criterion 2(b), the matching of expectation values of operators 
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• vertical or horizontal products 
■ diagonal products 
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Figure 4: Decay of type 2 couplings against length for k = 1 (left) and k = 
2 (right) (see table 1 for definitions) operators against length. Results are for 
3 RGTs at a w — 16 77 = .30. The solid line is the exponential fit, which is in 
agreement with Eq. 17. 



computed using different lattice sizes, we start by defining the quantitities 



a 



<o«> - (or 1 )) 



Q (l) = E«) 2 , (19) 

a 

which measure the quality of the matching for expectation values of opera- 
tors. A first insight is gained from the plot of the factor (see fig. 6). 
There is a clear minimum independent of aw for 77 ~ 0.3. This matching 
procedure appears to be a very precise criterion to fix the rescaling, but this 
a question to be considered later. What is more surprising from fig. 6 is 
the very good matching of observables for any value of aw including values 
for which we do not find evidence for a FP from criterion 2(a). However, 
criterion 2(b) is completely insensitive to the truncation of operators and to 
finite size effects, the systematic errors that affect criterion 2(a) the most, 
but is very sensitive to the canonical surface not being critical, which proba- 
bly accounts for the small systematic error in 77 (the exact value is rj = 0.25). 
A closer inspection, reveals that the quality of the matching is sensitive to 
aw- As apparent from fig. 7, a careful comparison shows that the approach 
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1 9 11 18 



0.25 


-0.6(2) 
-1.1(3) 


0.02(5) 
0.19(4) 


-0.03(7) 
0.15(8) 


0.002(20) 
-0.08(3) 


0.30 


-0.40(20) 
-0.87(21) 


-0.03(7) 
0.12(4) 


-0.03(7) 
0.14(8) 


-0.005(20) 
-0.09(2) 


0.40 


-0.40(20) 
-0.56(24) 


0.03(5) 
0.17(3) 


-0.03(6) 
0.12(7) 


-0.003(20) 
-0.06(2) 



Table 2: Example of couplings at 2 (up) and 3 RGTs(down) for aw = 40 that 
match poorly. 



1 18 21 30 



0.30 


0.67(7) 
0.60(6) 


0.014(6) 
0.00(1) 


0.001(5) 
0.015(6) 


0.0001(5) 
-0.0020(7) 


0.40 


0.52(6) 
0.38(6) 


0.012(5) 
-0.003(4) 


-0.001(4) 
0.012(4) 


-0.0000(4) 
-0.0015(6) 


0.50 


0.39(6) 
0.19(6) 


0.010(4) 
-0.003(3) 


0.001(4) 
0.010(4) 


-0.0000(4) 
-0.0010(4) 



Table 3: Example of couplings at 2 (up) and 3 RGTs(down) for a w = 8 with 
rather poor matching. 



to the FP is faster for a transformation at aw = 16. 

To cross-check criterion 2(b) further, we compared the matching of ob- 
servables on the next RGT, that is, observables computed in a 4 2 lattice. The 
number of operators available gets reduced, as some of them can no longer be 
accommodated in such small lattice, but for those that it exists, the matching 
should significantly improve as compared with the previous RGT. This is, 
indeed, corroborated from fig. 8, the matching in this case is close to perfect. 

5.3 The determination of the anomalous dimension 

Previous considerations have selected the value of aw = 16 as the one best 
compromising the 2 criterion of short-rangeness and fast approach to the 
FP. We now come to the determination of the rescaling factor, which as 
a byproduct allows us to compute the first critical exponent of the theory, 
the conformal anomalous dimension, see Eq. 15. We devised two different 
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Figure 5: Some chosen couplings at a w = 16 and i] = .25 for 2 and 3 RGTs. 



criterion to be meet. 

The first criterion 3(a) is the proper identification of a marginal eigenvalue 
in the T matrix. In our particular problem, this is complicated by the fact 
that in the cases were the two criteria we use to determine that a FP is 
reached, we clearly get two marginal directions. It seems natural to assume 
generally that the model has actually two marginal directions. Although we 
already ruled out some values of aw, we keep the discussion general for the 
time being and consider all values. The second eigenvalue of the T matrix 
as a function of rj is plotted in fig. 9 for 2 RGTs. Large values of aw are 
incompatible with the T matrix having a marginal eigenvalue, except for 
large values of the anomalous dimension, where the third eigenvalue is far 



18 




Figure 6: Dispersion in the matching of expectation values of operators, as 
compared from 3RGTs starting at 64 2 and 2 RGTs starting at 32 2 . 




Figure 7: Comparison of the quantity r a , as defined in Eq. 19 for aw = 8 
and aw = 16 as compared from 3RGTs starting at 64 2 and 2 RGTs starting 
at 32 2 . The RGT with a w = 16 produces a better matching. 
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Figure 8: Comparison of r a after 4/3(4 2 lattice) and 3/2 RGT (8 2 ) lattice 
at rj — 0.30. Only operators that can be accommodated in a 4 2 lattice are 
plotted. The matching is nearly perfect. 
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from one, and following previous discussions should be disregarded. This is 
reassuring, since for those values of aw we already know that RGT perform 
rather poorly. For smaller values of aw, we do get a one eigenvalue within 
statistical error bars in the range .10 < rj < .40. If a FP is reached, the 
marginal eigenvalue should be stable through the next RGTs, plotted in 
fig. 10. Also, for small values of rj, which criterion 3(b) will rule out as well, 
the third eigenvalue is clearly different from one. Within statitical error bars, 
the window of anomalous dimensions is further restricted to be in the range 
0.25 < rj < 0.30. 

The matching between operators, criterion 3(b) turns out to be the most 
sensitive criterion, a result already anticipated in fig. 6. For our favorite 
value aw = 16, all expectation values of operators match for rj = .30 within 
2a in the statistical error bars, for all 40 operators considered, see fig. 11. 
For comparison, we plot as well the relative error for values of the anomalous 
dimensions 77 = 0.10 and rj = 0.50 in which there is no matching at all. Let us 
recall that if we perform another RGT, the matching is perfect see previous 
fig. 8. This provides an important cross-check that we are in the vicinity of 
the IFP. Unfortunately, the final lattice (4 2 ) is too small and shows finite size 
effects in the critical exponents, as we will discuss. 

The matching criterion singles rj ~ 0.30 as the preferred value for the 
anomalous dimension, a result slightly off from the exact value, rj = .25. 
The matching at i] = .25 is also good but not as good as in 77 = .30. We 
will expand further on this issue in the discussion of the different sources of 
systematic errors. 

5.4 The critical exponents and the FP Hamiltonian 

The rescaling parameter has been computed and it corresponds to an anoma- 
lous dimension in the interval rj e (.25, .30), our results favoring values closer 
to rj = 0.30 than to rj = .25. 

In fig. 12 we plot the first eigenvalue as a function of the rescaling pa- 
rameter. The first eigenvalue is very sensitive to the value of the anomalous 
dimension. Statistical errors are relatively small, certainly very small when 
compared with the error bars in the the couplings. The first eigenvalue for 
the interval of anomalous dimensions considered is 

r] = 0.25 , Ai = 2.015(7) 
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Figure 9: The second eigenvalue of the T matrix as a function of the anoma- 
lous dimensions for different values of aw after 2 RGTs. 
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Figure 10: The second eigenvalue of the T matrix as a function of the anoma- 
lous dimensions for different values of aw after 3 RGTs. 
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Figure 11: Plot of the quantity r a as defined in Eq. 19 comparing the excellent 
matching at rj = 0.30 with no matching at all at r] = O.lO(top) and rj = 
.50(bottom). The final lattice is 8 2 . 
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Figure 12: First eigenvalue as a function of the anomalous dimension at 
aw = 16, final lattice is 8 2 . 
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Figure 13: Third eigenvalue as a function of the anomalous dimension at 
aw = 16, final lattice is 8 2 . 
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7] = 0.275 , Ai = 1.996(7) 
j] = 0.30 , Ai = 1.976(7) , 



(20) 



where the error is just statistical. Results are obtained after 3 RGTs, in a 8 2 
lattice, so that finite size effect errors are smaller than the statistical ones, 
as will be discussed. 

In fig. 13 the third eigenvalue is plotted and as already pointed out, it is 
compatible with one. 

We can compute the FP action. For consistency with our analysis, we 
will quote the results at rj = .30, although given the relative large error bars, 
it is essentially insensitive to r\ within the range (0.25,0.30). 



Operator 


K 


Operator 


K 


Operator 


K 


1 


0.16(12) 


8 


0.06(4) 


24 


-0.017(5) 


2 


0.99(8) 


9 


0.05(2) 


25 


-0.007(3) 


3 


-0.22(3) 


11 


0.07(4) 


27 


0.005(2) 


4 


0.03(5) 


12 


-0.06(2) 


33 


0.029(5) 


5 


-0.0016(3) 


15 


0.12(2) 


37 


-0.012(4) 


6 


-0.90(5) 


16 


0.03(1) 






7 


-0.17(4) 


21 


0.03(1) 







Table 4: FP Hamiltonian at a w = 16. Only the non-vanishing couplings are 
displayed. 



5.5 The different sources of errors 

We now analyze the different sources of errors that may appear. We treat 
each case separatedly. 

5.5.1 error e(l) 

This is the systematic error that may appear if some operators having sizeable 
coupling are not considered into Eq. 6. From previous results, we obtained 
that FP couplings decay at least exponentially, both with length and dis- 
tance, and as we parametrized the effective interactions generated along the 
flow within this assumption, we should not expect large systematic errors 
coming from truncation effects. 
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In the case of the evaluation of couplings, a typical plot of the evolution 
in the value of the coupling as more operators are considered is shown in 
fig. 14. Results are strongly sensitive to the inclusion of higher dimension 
operators but rather mildly to length. Furthermore, statistical error bars 
grow considerably with the inclusion of more operators. Eventually there is 
flat plateau, but only after enough number of operators are included. If we 
assume that the systematic error coming from the truncation, consistently 
with previous assumptins, decays exponentially, we estimate that is smaller 
than the statistical error. Although the truncation error is not completely 
negligible, statistical errors should be diminished to make it apparent. In 
other words, if larger statistics were available, higher dimensional operators 
should be considered. 

Concerning the evolution of eigenvalues, we do find a flat plateau which 
sets in with the inclusion of relatively few operators, see fig. 15 and fig. 16. 
For the case of eigenvalues we cannot find any evidence for a significant 
systematic error coming from the truncation. 

We conclude that the operators incuded are consistent within the statis- 
tics we considered. Errors coming from the truncation are smaller than the 
statistical one. 

5.5.2 error e(2) 

Finite size effects enter into the game when there are operators with sizeable 
coupling constants whose length is greater than the linear size of the lattice, 
L. Both the couplings and eigenvalues will then acquire a dependence on L. 
Of course this error should be more apparent as the lattice is smaller. Our 
results at a 4 2 do show finite size effects, but as our main conclusions are 
drawn from the analysis of a 8 2 lattice, we analyzed the effect in those. In 
fig. 17 we compare the couplings of the effective interaction generated after 2 
RGTs starting from a 32 2 and 64 2 lattices respectively. If there are finite size 
effects, these couplings should differ as they have different L dependences, 
irrespective of how we truncated the expansion (we truncate it in the same 
way in both lattices), and if we are reaching a FP( this is a question addressed 
by comparing successive RGTs). From fig. 17 we cannot detect a significant 
dependence on L within statistical errors. We conclude then, that finite size 
errors are within statistical ones if the smallest lattice used is a 8 2 one. 
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Figure 14: Evolution of coupling 6 as a function of operators included. 
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Figure 15: Evolution of the first eigevalue as a function of operators included. 
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Figure 16: Evolution of the second eigenvalue as a function of operators 
included. 
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Figure 17: Comparison of couplings at aw — 16, rj — .275 computed after 2 
RGTs starting from a 64 2 and 32 2 lattice respectively. 

5.5.3 error e(3) 

If the canonical surface is not exactly at criticality, the right infrared FP 
cannot be reached. The RG trajectories eventually flow towards the renor- 
malized trajectory instead, and the matching procedure should eventually 
deteriorate. Unfortunately, we cannot perform more RGTs to check if this is 
the case, and we do not have a more precise estimate of the critical couplings 
to compare with the results already obtained. From the matching criterion 
(precisely the method more sensitive to this error) rj = .30 is clearly preferred 
and the matching criterion is just marginally consistent with the exact result 
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r] = 0.25. This makes us suspect that there is a small systematic error of this 
kind in our calculations, which unfortunately we cannot estimate. 

5.5.4 error e(4) 

This error is associated with the transformation, and it is not peculiar to a 
finite lattice approximation. It appears when couplings in eigenvectors fail 
to converge. 

For the sake of clarity, in the appendix sect. B we present a simple model 
that can be solved exactly (the Gaussian model), where an explicit example 
of this problem is shown. 

An obvious criterion to detect this problem seems to compute the eigen- 
vectors, diagonalizing the T matrix, explicitly checking that they decrease 
as a function of length and dimension. For the first three eigenvalues (the 
2 and the two Is) that we already computed we do not find a decay of the 
eigenvector components as apparent as it was with the FP couplings. The 
next two eigenvalues at aw = 16 and r] G (0.25, 0.30) form a complex pair, 
so we think that those eigenvalues are affected by error e(4). As the eigen- 
vectors are complex as well, we do not know how to extract any consequence 
from the analysis of the coefficients. In any case, let us mention that finite 
lattice approximations to simpler models do also show complex pairs when 
the eigenvectors fail to converge, so it may be the case here as well. More 
statistics should clarify this issue. 

6 Conclusions 

6.1 Summary of the results 

In this paper we studied finite lattice approximations in models of uncon- 
strained spins. We discussed the properties that a renormalization group 
transformation should have in order that a finite lattice approximation may 
be organized into a controlled and systematic expansion, able to give accu- 
rate and rigorous results. We also discussed the difficulties that appear in 
computing the anomalous dimension of the field. 

As a non-trivial example, we studied the linear sigma model in two di- 
mensions in the domain of attraction of the Ising fixed point, using the Bell- 
Wilson RGT [4], which has a free parameter aw- The properties of the 
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RGT depend on aw and the goal was to properly identify the values of aw 
that optimally fullfill the two important properties relevant for us, the short- 
rangeness of the effective interactions generated and a rapid approach to the 
Ising fixed point. 

To address the issue of how short-ranged are the effective interactions 
essential use has been made of the Schwinger-Dyson equations Eq. 7 criterion 
1(a). Overall the method performs well, but the statistical error bars are 
relatively large, which in some cases limits the accuracy of the results. In 
any case, we have provided convincing evidence that the decay ansatz Eq. 12 
is well satisfied (see fig. 2, fig. 3 and fig. 4), which in turn provides the 
justification for the expansion we propose for the effective interactions. Our 
results clearly show that the smaller the value of aw the better the decay 
properties and that decay with dimension is slower than with distance. If 
more statistics would be available, the expansion should include more higher 
dimensional operators than higher length. 

Although the smaller aw, the more short-ranged the fixed point action 
is, we found that for slightly higher values (aw ~ 16) the canonical surface 
is closer to the corresponding FP. To reach this conclusion we computed the 
couplings generated after successive RGTs and checked whether they agree, 
criterion 2(a). This was helpful to rule out high values of aw, and to provide 
strong evidence for a FP at aw = 16(77 = 0.25) (see fig. 5). We also used the 
matching criterion of observables 2(b) and found only a tiny dependence on 
aw, enough to show that aw = 16 is the optimal case (see fig. 7). 

Having identified the value aw = 16 as the best suited for a finite lattice 
approximation, we computed the rescaling of the field, and found an anoma- 
lous dimension i] G (0.25,0.30). Although compatible with the exact result 
i] = 0.25, our results favor values closer to rj = .30. The first criterion used 
to address this issue, the eigenvalue being one, led to a window of acceptable 
anomalous dimension (0.25, 0.30) after a further assumption that there are 2 
marginal directions. This further assumption comes as a particularity of the 
two-dimensional linear sigma model and should not show up in other models. 
The matching criterion 3(b), which turned out to be a very sensitive criterion 
(see fig. 11), clearly favored rj ~ 0.30. 

After all this previous study, we present the final results for the critical 
exponents, 

77 E (0.25,0.30) (21) 



30 



v = 1.00(2) 
A 2 = 1.0(4) 
A 3 = 1.0(4) , 



(22) 
(23) 
(24) 



where the error merely indicates the uncertainty in the determination of the 
anomalous dimension, which in turn, depends on the statistical error and 
the lack of criticality, as explained. We would like to emphasize that all the 
rest of systematic errors have been considered and shown to be smaller (see 
the discussion at the end of the last section). We also computed the FP 
Hamiltonian and found that it may be parametrized with a few operators 
(see table 4). Finally, we provide arguments that within a Bell- Wilson RGT 
one can only compute a finite number of exponents, the rest are not accessible 
to numerical methods due to error e(4). Most likely, this number reduces the 
accessible eigenvalues to the ones we quote, although it is possible that this 
is just a reflection of the modest statistics used. 

6.2 Perspective and outlook 

By optimizing the transformation used we have been able to get a great deal 
of information, including precise calculations for critical exponents, FPs, etc. 
with modest statistics, and more importantly, to have under control the dif- 
ferent source of errors. There are, however, some problems that we should 
emphasize. Concerning the determination of the anomalous dimension, our 
results favor values of rj ~ .30, although rj = 0.25 is still acceptable. Higher 
statistics and a more accurate determination of the critical surface should 
resolve this issue. Another important point concerns real space renormal- 
ization group methods as they seem to present a limitation concerning the 
number of, generally irrelevant, critical exponents accessible within finite lat- 
tice approximations. Deeper understanding of the properties of RGTs are 
necessary with the ultimate goal of proposing RGTs with better convergence 
properties and good short-rangeness behaviour, as is discussed in the final 
appendix. 

There are similar MCRG calculations for models of unconstrained spins 
using momentum RGTs [16]. As discussed, RGTs in momentum space have 
different properties and are not directly comparable. Closer works to com- 
pare ours with are MCRG calculation for 3d Ising models using real space 
RGTs [17], [18] and the most recent update [19]. Although the large statistics 
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and big lattices reported in those papers are completely unattainable with 
our resources, at a more conceptual level there are interesting analogies we 
can draw. The RGT transformation used in those works is a majority rule, 
which is the spin-constrained version of the Bell- Wilson RGT at a w = oo. 
We have seen that this is a very unfavorable case concerning the short-ranged 
properties of the effective interactions and the approach to the FP, so it may 
be the case as well for the constrained model (which can be obtained as the 
A — > oo case of the unconstrained one). Let us recall that this is not at 
odds with the good matching for expectation values of operators. We find 
a surprisingly accurate matching for 3 RGTS even at aw = oo, despite that 
other criteria shows that we do not reach a FP after 3RGTs. Our best value 
for the matching singles out an anomalous dimension slightly off from its 
exact value. Similarly, in 3d Ising MCRG calculations the anomalous dimen- 
sion comes slightly off from its more accepted value [20]. Secondly, only two 
even parity eigenvalues may be computed, a result which seems in agreement 
with the difficulties mentioned concerning the convergence of higher irrel- 
evant eigenvalues. It would be interesting to perform similar calculations 
using more elaborate RGTs. 

Although, as discussed, some questions demand further clarification, we 
have shown that with a suitable knowledge of the RGT, finite lattice approx- 
imations provide a systematic expansion to perform calculations in the RG. 
We hope that this paper will help to arouse and inspire further work in this 
field, which will allow application of these methods to more general models. 
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A The Strong coupling expansion 



We determine the CS by means of a high temperature expansion (HTE) of 
the susceptibility, 

oo 

X = EX? ) «S ( 25 ) 

i=0 

where depend on A through integrals like J ra (A) = d(fxf) n e~ s ^ . Baker 
and Kincaid [21], have tabulated its first 11 coefficients, and in such a case, 
only even integrals up to n = 12 have to be considered. 

Had we computed the whole series, the critical coupling is obtained from 
the ratio method [22] , 

K <=M = J™/"' r n = Xi n_1) /Xi n) > ( 26 ) 

but in our case must be extrapolated from finite n ratios. In such a case 
corrections to the preceeding formulae are correlated with the form of the 
closest singularity of the susceptibility to the origin in the complex k plane. 
It is assumed that, 

^MA)(l + ^ + 0), (27) 

plus some exponential corrections exp(— na) arising from singularities farther, 
provided that within the disk \k\ = n c (\) there are no further singularities. 
This is, unfortunately not the case for a PSQ (plane square) lattice (as it is 
always for loose-packed lattices), due to the presence of an antiferromagnetic 
singularity at k — —k c (X), which brings in nonanalytical oscillating power 
law corrections in Eq. 27, of the type, 

~ ^ (28) 

$ is the critical exponent associated to the antiferromagnetic singularity. 

The ratios in Eq. 26, of the expansion in Eq. 25, are nearly constant with 
a monotonic decreasing oscillating trend. Using a double ratio method [22], 
a rough estimate of the critical coupling may be obtained K re . We define 
then a new variable through, 

z = rfV. (29) 

Kre 
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so that we send the antiferro magnetic singularity almost to oo, and thus the 
susceptibility expressed in this new variables eliminates the oscillating term 
in Eq. 28 as it picks up a large exponential correction. 

We can extrapolate now by using Eq. 27 for n > 5 and we estimate 
the error from comparing extrapolations including different invers powers in 
n. We also analyzed the ratio method if insted of K re , we input the new 
better estimate, but nothing is gained within error bars. In any case, Just 
to compare with known results, we obtained for A = oo (Ising limit), 

K c = 0.4407(2), (30) 

in agreement with the exact result K c = 0.4406868. 



B The Non-convergence of eigenoperators 

In this appendix we show in an explicit example how exact eigenoperators of 
the T-matrix may not converge. We explictly show as an example the Gaus- 
sian Model, defined as the most general translational-invariant quadratic 
Hamiltonian 

^ = ^EEP«0(n)0(n + r), (31) 

n r 

where the sum is over sites of the corresponding lattice. Under /c-iterations 
of a linear RGT the model transforms into itself, 

W* = ^EEP*( r )#(«)#(« + r) • (32) 

" n r 

It is possible to derive an exact RG-equation for the RGT used in this paper 
[4] , which is not difficult to generalize to a completely general linear RGT [6] . 
We can actually solve the equations and find the exact FP and eigenoperators. 
We just quote the final result for the j-eigenoperator [6] 

+oo -y d 

&jP(p) = P*(p) 2 9j(p) , 9j(p) = E ( p + 27r/)( 4 - 2 i) fl^Ov + 27r W > ( 33 ) 

where the equation is expressed in momentum space, j = 0, 1, • • •, the di- 
mension of the system is d (d — 2 in our case), and p* is the FP value of 
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the couplings which is not directly relevant for our problem and will appear 
elsewhere. The function g(x) depends on the RGT. For the Bell- Wilson RGT 
Eq. 13 we have [4] 



g(x) = 



sin^f) 
(f) 2 



(34) 



In this case, eigenoperators Eq. 33 are well defined for j = (relevant), j = 1 
(marginal), j = 2 (irrelevant), but the sum in Eq. 33 fails to converge for j > 
2. Any numerical method will necessary fail to get this eigenoperator right, 
as it has just a formal meaning. This conforms very well to the numerical 
results reported in [4]. 

To remedy this situation, one has to consider other RGTs, which have 
overlapping cells (see fig. 18). Let us mention that these transformations 
have appeared in the context of perfect actions [5]. For such a case one gets 
a g-function 

sm 4 (x) 



g{x) = 



(f) 4 



(35) 



so now, for example, from Eq. 33 the j = 3 eigenoperator does converge, 
opposite to what it happens with Bell-Wilson RGT. In other words, the 
eigenvalue with j = 3 becomes now accessible to numerical methods. In fact, 
the RGT in fig. 18 is just the simplest case of an overlapping RGT. One may 
define more elaborated RGTs so that the g— function Eq. 34 becomes more 
convergent [6]. 



C=2 



Figure 18: The simplest case of overlapping RGT. 
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Let us recall that RGTs in momentum space are an extreme case of over- 
lapping RGTs. In such a case, it is easy to prove that eigenoperators are 
perfectly convergent for any value of j. However, this is a very unfavorable 
case concerning the short-rangeness of the effective actions generated, be- 
cause the decay ansatz Eq. 12 gets replaced by an algebraic decay It seems 
then, that there is some kind of uncertainty principle concerning RGTs; One 
may have very well behaved RGTs concerning how short-ranged the effec- 
tive interactions are at the expense of some eigenoperators having just a 
formal meaning, or one may choose RGTs having very well defined set of 
eigenoperators at the expense of having poor short-ranged behaviour of the 
Hamiltonian generated along the RG flow. 
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